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PROBLEMS  ON  THE  GEOGRAPHICAL  DISTRIBUTION  OF  THE 
ATMOSPHERIC  ELECTRIC  FIELD  INDUCED  BY 
THUNDERSTORM  ACTIVITIES* 

by 

Lee  Chi-chen  and  Lin  Chi -zhang 

Atmospheric  electricity  has  been  a  formal  subject  of  study  for 
more  than  two  hundred  years.  The  wealth  of  observational  evidence 
accumulated  during  the  long  period  in  the  past  reveals  that  the  electri¬ 
cal  phenomena  observed  on  fine  days  are  distinctly  different  from  those 
associated  with  bad  weather.  The  field  of  electric  potential  in  regions 
of  fine  weather  varies  uniformly  with  the  activities  of  thunderstorms  on 
a  global  scale.  Electric  phenomena  in  the  atmosphere  over  good  and 
bad  weather  regions  are  thus  interrelated  [1], 

The  world-wide  nature  of  atmospheric  electricity  is  generally  well 
known.  However,  a  satisfactory  theoretical  explanation  on  the  cause  of 
electric  activities  in  the  atmosphere  is  still  lacking. 

Two  relatively  important  theories  have  been  proposed,  namely, 

Wilson's  theory  of  spherical  condenser  and  the  theory  of  Frenkel.  Wilson 

[2]  has  considered  that  the  earth  and  the  ionosphere  form  two  poles  of  a 

spherical  condenser  and  that  thunderstorms  are  equivalent  to  a  dynamo. 

Due  to  over  simplification  of  the  problem,  his  treatment  is  unable  to 

account  for  the  observed  fact  that  the  distribution  of  atmospheric  charges 

occur  in  the  troposphere  instead  of  in  the  ionosphere.  Moreover,  his 

model  is  too  crude  for  a  rigorous  analysis  to  be  made  by  means  of  mathematical 
* 
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physics.  In  Frenkel's  theory  [1,  3],  a  zero  potential  at  the  earth's 
surface  is  assumed  or  the  total  electric  charge  or  the  entire  surface 
of  the  earth  is  taken  to  be  zero.  Consequently,  Hs  solution  to  the 
problem  of  an  infinitely  large  plane  conductor  earthed  to  zero  electric 
potential,  or  of  a  neutral  conductor  deviates  markedly  from  actual 
atmospheric  conditions  and  the  results  cannot  be  regarded  as  conclu¬ 
sive. 

In  t^fs  paper,  the  electromagnetic  equation  in  spherical  coordinates 
is  , first  derived  with  due  consideration  on  the  actual  electrical  conduc¬ 
tivity  of  .the  atmosphere  and  the  correct  selection  of  conditions  for  solu¬ 
tion.  The  distribution  of  the  atmospheric  electric  field  induced  by  thun¬ 
derstorm  activities  is  then  evaluate d.^The  results  are  found  to  be  better 
than  those  given  by  Wilson  and  Frenke^L 


I.  THE  BASIC  EQUATION  AND  CONDITIONS  FOR  SOLUTION 

For  the  analysis  of  the  electric  processes  in  the  atmosphere  the 
Maxwell  electromagnetic  equations  may  be  expressed  in  conventional 


symbols  as: 


V  •  E  *«  4* p, 

V  X  E  =  -  - 

i 

V  •  B  =■  0, 


J_  9B 
c  di  5 


vxB=-j+i 

c  c  dt  ) 

The  conservation  of  electric  quantity  requires  that 

If  "  “V  '  i*  (2) 

In  addition  to  the  conduction  current  j  which  obeys  Ohm's  law, 
a  current  j  may  also  be  induced  by  an  external  electromotive  potential 
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in  the  atmosphere.  The  generation  and  separatum  of  electric  charge 
in  cloud  give  rise  to  electric  currents  which  do  not  obey  Ohm's  law. 
Therefore,  we  have 


h  +  J  , 


(3) 


h  -  AE,  {4) 

where  \  denotes  the  atmospheric  conductivity. 

In  the  application  of  the  above  equations  to  the  analysis  of  elec¬ 
tric  processes  in  the  atmosphere,  it  is  more  convenient  to  use  a 
spherical  polar  coordinate  system  (r,  0,  <p)  with  the  origin  at  the  center 
of  the  earth.  Calculations  are  made  on  the  following  three  basic  assump¬ 
tions: 

(a)  Quasi -stability:  As  the  relaxation  time  of  the  electrical  pro¬ 
cesses  in  the  earth-atmosphere  system  is  relatively  small  [4],  it  is 
considered  that  a  quasi -equilibrium  state  may  exist  at  any  particular 
time;  thus  we  have  d/dt  =  0. 

(b)  According  to  observations,  the  electrical  conductivity  varies 
exponentially  with  height  [2,  5]  and  is  given  by 

A(r,  0,  ,P)  -  U0,  ric-W-'J,  (5) 

where  r0  denotes  the  radius  of  the  earth. 

(c)  Horizontal  quasi-uniformity:  Electrical  conductivity  in  the 
atmosphere  is  essentially  a  function  of  height.  Although  its  value  may 
vary  from  place  to  place,  the  horizontal  gradient  is  much  smaller  than 


i 

( 

i 


the  vertical  gradient.  Thus  terms  with  the  horizontal  gradient  of 
X  in*  the  equation  may  be  neglected. 

The  first  assumption  gives  dB/dt  *  0  and  vxE  =  0.  Hence  the 
electric  potential  V  may  be  introduced  as  E  «  -  v  V.  Since  dp/dt  =  0, 
we  have 

V  -  j  -  -XV2V  -  VX  •  W  +  V  •  j'  -  0.  \ 


The  generation  and  separation  of  electric  charge  in  cloud  form  one 

type  of  "external1’  electromotive  potential  and  give  rise  to  an  electric 

current  j',  which  maintains  several  centers  of  electric  charge  in  cloud. 

A  positive  (negative)  charge  center  Q.,  at  (r^,  0..,  <p^)  acts  as  a  sink 

(source)  of  j'  but  as  a  source  (sink)  of  j  .  This  leads  to 

c 


-  V  •  3  =  V  •  }'  -  -4 7t0,l,o(r  -  r,)8(0  -  0,)$(<p  -  <p,), 


(7) 


where  Q..  and  X^  denote  the  quantity  of  electric  charge  and  the  electric 
conductivity  at  a  particular  charge  center  respectively,  while  6(x)  repre¬ 
sents  a  unit  pulse  function.  At  this  stage,  we  may  treat  each  charge 
center  as  a  point  charge  and  then  extend  the  treatment  to  a  more  general 
distribution  of  electric  charge. 

By  substituting  Equation  (7)  into  (6)  and  considering  the  second  and 
third  assumptions,  we  obtain 

V  ^  +  a  ~dr  *”  ""  **Qig(r  —  ri)s(.d  —  0,)s((p  —  <p,),  {g) 

where  ^  represents  a  summation  of  all  electric  charge  at  the  individual 

r.hnr<r<>.  i  in  thfi  af.rnnRaherft. 
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Since  the  earth  is  a  good  conductor,  the  earth's  surface  r  -  v0 
may  be  taken  as  an  equipotential  surface.  If  we  assume  as  usual 
that  the  electric  potential  is  zero  at  infinity,  then  the  potential  at  the 
surface  of  the  earth  V0  remains  to  be  determined,  Frenkel  [3]  has 
taken  V0  as  zero  and  this  is  one  of  the  factors  which  make  his  approach 
unsuccessful.  According  to  the  first  assumption,  the  time -variation  of 
the  total  electric  charge  Qeartk  at  the  earth's  surface  may  be  expressed 
as  ®^earth^®*  *  Thus  sum  of  all  electric  currents  at  the  earth's 
surface  S^,  should  be  equal  to  zero. 

The  above  formulation  represents  an  attempt  of  finding  a  solution 
to  the  following  boundary-value  problem: 


V>V  +  “  fr  “  “  **QMr  -  r,)e('p  -  >p,)8(,0  -  0,)v 


V( 0,  i /<)  “  F,,  iim  J7(r,  0,  ij>)  «=»  0, 

r  »*• 

I'-f 


(8) 

(9) 

(10) 


n.  SOLUTION  OF  EQUATION  AND  CONCLUSION 
By  making  use  of  the  second  assumption  and  the  fact  that  a/r  «  az/4, 

a(6.  <?)/r  r  \ 

we  may,  after  a  change  of  variables  V(r,  0,  <p)  *  u(r,  0,  <p)e  2  '  0  , 

transform  Equation  (8)  into  a  typical  elliptic  partial  differential  equation 
which  gives  a  solution  to  the  above  boundary-value  problem  [  5]  in  units 
of  coulomb -volt-km  as: 


\r,  V,  <p)  10’  2_j  i— 

<  h 


;  ! 


i  : 


1  1 
5 

i  ' 


I"  - 

^0  ^  1  I  t»  —  \  fa«  t 

1.  R, 

,  +  »V  (11) 

rl  Kf  J 

'V 

0 
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where 

r)  ==»  — ,  i?/  “  \f  r1  +  r)  —  2rr/cesy,  R\  *=»  r1  4-  r}'  —  2rrJ  cos  y 
ri 

cosy  *=*  cos 0  cos  +  sin  6  sin  0/  •  cos  (qo  —  rp/) 

while  R.  and  R’.  denote  the  distance  from  (r.,  9.,  9).)  and  (r1.,  9^,  <p  to 
the  point  of  observation  (r,  9,  <p )  respectively. 

Direct  verification  reveals  that  Equation  (11)  satisfies  Equation  (8). 
The  boundary  condition  of  Equation  (9)  is  satisfied  when  rl/jq,  .1,  because 
the ifirst  term  in  the  above  expression  is  approximately  zero  at  r  =  r0  and 
the  second  term  in  Equation  (11)  becomes  V0  at  r  =  r0 . 

The  first  term  on  the  right-hand  side  of  Equation  (11)  represents  the 
electric  field  induced  by  the  charge -centers  in  cloud,  while  the  second 
term  denotes  the  charge  distribution  of  the  conducting  sphere  (the  earth). 
Thus  they  may  be  denoted  by  Vi(r,  0,  (p)  and  V2  (r,  0,  )  respectively.  In 

order  to  determine  the  coefficient  V0  in  V2(r,  0,  <p),  we  may  substitute 
Equation  (11)  into  Equations  (10)  to  get 


0,  +'t>  ,r 

—  — - - - = 

0,  ip)n(d,  (p)iiS 


where  Rt  denotes  the  total  resistance  of  the  earth-atmosphere  system  and 
Ij.  the  sum  of  the  conduction  currents  of  all  bad  weather  regions  over  the 
world,  currents  induced  by  lightning  and  those  released  by  sharp  points 
and  precipitation.  According  to  References  [2]  and  [6],  1^  is  about  1530- 
1800  amperes  and  R^  about  200  ohms,  giving  V0  equal  to  -3.06  to  -3.6  x 
X  10s  volts. 


The  above  results  lead  to  the  following  points  for  discussion: 

(a)  Differentiation  of  V(r,  0,  <p)  with  reapect  to  r  gives  the  poten 


tial  gradient  at  the  surface  of  the  earth  as: 
dV(r,  0,  ip)\  _  ,  v>  „  / _ \ _ 


r,T 

hJL\l 

1  /?) 

\2 

h  J  J 

The  results  of  the  computed  values  of  V{r,  0,  <p)  and  9V/9r  _  are 

r  -  i'0 

shown  in  Figure  1. 


Distance  (km) 
(*) 


<b) 

Figure  1 


(a)  Distribution  of  electric  potential  in  the 
region  of  a  thunderstorm 
(Unit:  104  volts.) 

(b<  Distribution  of  the  surface  electric  field 
in  the  regie::  of  a  thunderstorm 
(Unit:  volts/cm.) 
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(b)  The  fu*3t  terms  in  Equations  (U)  and  (13)  represent  the  field 
of-electric  <charge  in  cloud.  It  is  noted  that  one  additional  factor 

exp  |+  (r.  -  r0 ) - is  involved  in  comparison  with  similar 

terms  for  a  static  environment  problem.  Thus  it  is  shown  that  the 
effect  of  a  point  source  decreases  with  distance  more  rapidly  in  an 
atmosphere  with  non-uniform  conductivity  than  in  a  non-conducting 
atmosphere.  This  point  is  worthy  of  attention  because  its  negligence 
may  lead  to  an  underestimation  of  the  values  of  Q.  in  assessing  the 
potential  gradient  at  the  earth's  surface  in  the  vicinity  of  thunderstorms. 

(c)  The  distribution  of  electric  potential  in  fine  weather  regions 

at  a  distance  of  more  than  10  km  from  a  point  source  may  be  determined 
by  the  second  term  V2(r,  9,  <p )  in  Equations  (9)  and  (11),  which  repre¬ 
sents  the  field  of  negative  charge  of  the  earth.  It  may  be  noted  that 
V2  0  when  r  -*■  »  and  V  -*•  V0  <0  when  r  -*■  0,  indicating  that  equi- 
potential  conditions  exist  at  the  earth's  surface  and  at  the  upper  levels. 

The  potential  difference  is  given  by  AV  =  -  V0  -  I^R^.,  i.  e.  >  the  upper-air 
electric  potential  is  higher  than  that  at  the  earth's  surface  by  300  to  400 
kilovolts.  At  the  earth's  surface,  we  have  9V2/9r  =  -  aV0  ~  108  volts/m. 
These  values  are  found  to  be  in  a  good  agreement  with  observations. 

(d)  The  stronger  the  global  thunderstorm  activity  the  larger  is  the 

value  of  It,  giving  a  higher  value  of  -V0  correspondingly.  This  gives 
rise  to  the  so-called  "uniform  variation"  in  the  electric  field.  Since  the 
value  of  a(9,  (p  )  may  vary  with  space  and  time,  changes  in  the  potential 
gradient  at  the  earth's  surface,  9V2/9r  _  =  a(9  ,  <p)V0  are  only  uniform 

over  the  sea  where  a  is  approximately  constant. 


1 

f 


(e)  In  fine  weather  regions,  we  have  V  =  r-r°  ^  for 

the  electric  potential,  -  -aV0  e  ^  r°  ^  >  0  for  the  potential  gra¬ 
dient,  p  =  “•  e  a(r~r°  ^  >  0  for  the  charge  density  and 

Qearth  =  f  dS  <  0  for  the  total  char8e  at  the  earth's 

surface.  Thus,  the  electric  potential  increases  exponentially  with 
height,  while  the  potential  gradient  and  the  charge  density  decrease 
exponentially  with  height.  These  are  in  good  agreement  with  average 
observed  conditions.  In  addition,  the  results  of  computations  show 
that  the  distribution  of  positive  charge  occurs  in  the  troposphere  with 
negative  charge  near  the  earth's  surface. 

In  the  above  calculations,  the  charge -centers  in  cloud  are  treated 
as  point  sources.  The  derived  results  may  easily  be  generalized  to 
cover  a  continuous  distribution  of  electric  charge  in  the  form  of 
y(.r,  0,  tp)  =  V'e-WX'-r.)  4- 

(i4) 


where  r  denotes  the  position  vector,  while  p(ri)  and  t.  represent  the 

til 

charge  density  distribution  of  the  i  charge -center  and  the  total  volume 
respectively.  In  discussing  the  electric  field  in  the  vicinity  of  a  charge- 
center,  it  is  necessary  to  determine  accurately  the  distribution  of  charge 
density  p{r/}  [  7].  However,  the  requirement  does  not  affect  the  fore- 
going  conclusive  remarks. 


0 
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